Andrew算法(我确实不懂Graham) 您所在的位置:网站首页 andrew s Andrew算法(我确实不懂Graham)

Andrew算法(我确实不懂Graham)

#Andrew算法(我确实不懂Graham)| 来源: 网络整理| 查看: 265

先解释一下:这两个算法分别都是凸包问题的算法,然后Andrew是Graham的变种,速度更快,更稳定,非常优秀,介于我已经把Graham写的莫名其妙的WA了,所以我选择了这种算法!

我认为在这里,还是有必要向大家在这里先普及一下什么是凸包:

计算几何凸包

凸包:给你n个散落的点,让你求出最小的凸多边形将所有的点包括起来,或者点在边上。用到的算法是Graham或Andrew!

这里给一个链接: Graham算法详解

下面给出Andrew算法的思路 First

首先先按照x坐标大小或y坐标大小进行排序(如果x坐标一样,y坐标就从小到大排序或如果y坐标一样那么x坐标就从小到大排序)

Second

然后进入程序的主干部分,先说一下Andrew主干的大体思路,我们分两次来求这个凸包,先从左到右一遍,再从右到左一遍(或先从下到上一遍,再从上到下一遍)首先我们一定要明白第n-1个点一定会在第一遍时进入凸包栈内(看了上面链接的朋友都应该知道这个栈是如何操作的,这里不再赘述)(因为n个点是从0n-1),所以第二遍的时候不必从n-1开始,从n-20开始就可以了(代码上会有体现)然后就完了!

现在我们来详细讲一下如何实现Second的操作

我们要实现找凸包,那么就必须找到最外层的点,这里就要使用叉积进行判断(向量a叉向量b=a.x×b.y-b.x×a.y)如果为正a在b的右边反之在左边(题目中因为我们只能定义一个基准坐标系,所以为了实现这个功能我们就必须找参照点,参照点作为临时原点)代码中xmult会体现。

然后就差不多了!

下面就是代码了:

#include #include #include #include using namespace std; struct point{ int x,y; }; bool cmp(point a,point b) { if(a.y==b.y&&a.xbasic&&xcross(node[i],node[num[top]],node[num[top-1]])) top--; top++; num[top]=i; } double s; s=0.0; for(int i=1;i


【本文地址】

公司简介

联系我们

今日新闻

    推荐新闻

    专题文章
      CopyRight 2018-2019 实验室设备网 版权所有